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Abstract 

Simple SUSY GUT models based on the gauge group SO (10) require t — b — r Yukawa 
coupling unification, in addition to gauge coupling and matter unification. The Yukawa 
coupling unification places strong constraints on the expected superparticle mass spec- 
trum, with scalar masses ~ 10 TeV while gaugino masses are quite light. A problem 
generic to all supergravity models comes from overproduction of gravitinos in the early 
universe: if gravitinos are unstable, then their late decays may destroy the predictions 
of Big Bang nucleosynthesis. We present a Yukawa-unified 50(10) SUSY GUT scenario 
which avoids the gravitino problem, gives rise to the correct matter-antimatter asymmetry 
via non-thermal leptogenesis, and is consistent with the WMAP-measured abundance of 
cold dark matter due to the presence of an axino LSP. To maintain a consistent cosmology 
for Yukawa-unified SUSY models, we require a re-heat temperature Tr ~ 10 6 — 10 GeV, 
an axino mass around ~ 0.1 — 10 MeV, and a PQ breaking scale f a ~ 10 12 GeV. 
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1 50(10) SUSY GUTs and Yukawa unification 



Grand unified theories (GUTs) are amongst the most compelling ideas in theoretical physics. 
Their beauty is only enhanced via a marriage to supersymmetry (SUSY). The SU(5) theory [I] 
unifies the Standard Model (SM) gauge symmetries into single Lie group, while explain- 
ing the ad-hoc hypercharge assignments of the SM fermions, and successfully predicting the 
m,b/m T ratio. Adding SUSY to the SU(5) theory stabilizes the hierarchy of interactions, but 
also receives experimental support from the celebrated unification of gauge couplings at scale 
M GUT ~ 2 x 10 16 GeV. 

The SO (10) SUSY GUT theory has even further successes[2]. For one, it explains the ad-hoc 
anomaly cancellation within the SM and SU(5) theories. Further, it unifies all matter of a single 
generation into the 16-dimensional spinor representation ^(16), provided one adds to the set of 
supermultiplets a SM gauge singlet superfield Nf (i = 1 — 3 is a generation index) containing 
a right-handed neutrinoQ Upon breaking of SO (10), a superpotential term / 3 ^M^N^Nf is 
induced which allows for a Majorana neutrino mass which is necessary for implementing the 
see-saw mechanism for neutrino masses [1]. In addition, the 5*0(10) theory allows for unification 
of Yukawa couplings of each generation. This applies calculationally especially to the third 
generation, where in simple SO (10) SUSY GUTs, we may expect t — b — r Yukawa coupling 
unification in addition to gauge coupling unification at scale Q = Mqut® E]. 

In spite of these impressive successes, GUTs and also SUSY GUTs have been beset with a 
variety of problems, most of them arising from implementing GUT gauge symmetry breaking 
via large, unwieldy Higgs representations. Happily, in recent years physicists have learned 
that GUT theories- as formulated in spacetime dimensions greater than four- can use extra- 
dimension compactification to break the GUT symmetry instead[7]. This is much in the spirit 
of string theory, where anyway one must pass from a 10 or 11 dimensional theory to a 4-d 
theory via some sort of compactification. 

Regarding Yukawa coupling unification in 50(10), the calculation begins with stipulating 
the b and r running masses at scale Q = Mz (for two-loop running, we adopt the DR reg- 
ularization scheme) and the t-quark running mass at scale Q = m t . The Yukawa couplings 
are evolved to scale Q = Msusy, where threshold corrections must be implemented [8], as one 
passes from the SM effective theory to the Minimal Supersymmetric Standard Model (MSSM) 
effective theory. From Msusy on to Mqut-, Yukawa coupling evolution is performed using 
two-loop MSSM RGEs. Thus, Yukawa coupling unification ends up depending on the complete 
SUSY mass spectrum via the t, b and r self-energy corrections. 

In this letter, we adopt the Isajet 7.75 program for calculation of the SUSY mass spectrum 
and mixings [9] and IsaReD[10] for the neutralino relic density. Isajet uses full two-loop RG 
running for all gauge and Yukawa couplings and soft SUSY breaking (SSB) terms. In running 
from Mgut down to M wea k, the RG-improved 1-loop effective potential is minimized at an 
optimized scale choice Q = y/mi L mi R) which accounts for leading two-loop terms. Once a 
tree-level SUSY/Higgs spectrum is calculated, the complete 1-loop corrections are calculated 
for all SUSY/Higgs particle masses. Since the SUSY spectrum isn't known at the beginning of 
the calculation, an iterative approach must be implemented, which stops when an appropriate 
convergence criterion is satisfied. 

1 Hcre, we adopt the superfield "hat" notation as presented in Ref. [3]. 
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Yukawa coupling unification has been examined in a number of previous papers [SJ EJ EH [E2J 
IT3"| [H] . The parameter space to be considered is given by 

mi6, raw, M 2 Dl mi/2, A , tan/3, sign(n) (1) 

along with the top quark mass, which we take to be m t = 171 GeV. Here, m\Q is the common 
mass of all matter scalars at Mqut, mio is the common Higgs soft mass at Mqut and Mp 
parameterizes either D-term splitting (DT) or Higgs-only soft mass splitting (HS). The latter 
is given by m 2 H = mf =F 2M|>. As in the minimal supergravity (mSUGRA) model, min 
is a common GUT scale gaugino mass, A Q is a common GUT scale trilinear soft term, and 
the bilinear SSB term B has been traded for the weak scale value of tan/3 via the EWSB 
minimization conditions. The latter also determine the magnitude (but not the sign) of the 
superpotential Higgs mass term /i. 

What has been learned is that t — b — r Yukawa coupling unification does occur in the MSSM 
for fi > (as preferred by the (g — 2) M anomaly), but only if certain conditions are satisfied. 

• The scalar mass parameter m 16 should be very heavy: in the range 5-20 TeV. 

• The gaugino mass parameter m^j should be as small as possible. 

• The SSB terms should be related as A% = 2mf Q = ^m\ G , with A = — 2mi 6 (in our sign 
convention). This combination was found to yield a radiatively induced inverted scalar 
mass hierarchy (IMH) by Bagger et al. [15] for MSSM+right hand neutrino (RHN) models 
with Yukawa coupling unification. 

• tan (3 ~ 50. 

• EWSB can be reconciled with Yukawa unification only if the Higgs SSB masses are split 
at Mqut such that rn 2 Hu < m 2 Hd . The HS prescription ends up working better than DT 
splitting [I3JE2]. 

In the case where the above conditions are satisfied, then Yukawa coupling unification to 
within a few percent can be achieved. The resulting sparticle mass spectrum has some notable 
features. 

• First and second generation matter scalars have masses of order m 16 ~ 5 — 20 TeV. 

• Third generation scalars, and fi are suppressed relative to tuiq by the IMH mechanism: 
they have masses on the 1 — 2 TeV scale. This reduces the amount of fine-tuning one 
might otherwise expect in such models. 

• Gaugino masses are quite light, with nig ~ 350 — 500 GeV, ~ 50 — 80 GeV and 
m lT/ - 100 - 150 GeV. 

The sparticle mass spectra from S'O(IO) SUSY GUTs shares some features with spectra 
generated in "large cutoff supergravity" or LCSUGRA, investigated in Ref. [TB]. LCSUGRA 
also has high mass scalars- typically with mass around 5 TeV- and low mass gauginos. The 
.50(10) SUSY GUT models are different from LCSUGRA in that they have a large A , with 
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A ~ — 2m 16 , and a \i term of around 1-2 TeV. This means 5*0(10) SUSY GUTs have a 
dominantly bino-like Z\ state, whereas the LCSUGRA authors adopt the mSUGRA model 
focus point region, which has a mixed higgsino-bino Z\ state. The latter can easily give the 
measured abundance of cold dark matter (CDM) in the form of lightest neutralinos. 

Since the lightest neutralino of 5*0(10) SUSY GUTs is nearly a pure bino state, it turns 
out the neutralino relic density Q^h 2 is calculated to be extremely high, of order 10 2 — 10 4 . 
This conflicts with the WMAP-measured value[T7], which gives 

tt CD Mh 2 = P CDM/Pc =0.11ltHi5 (2a). (2) 

where h = 0.74 ± 0.03 is the scaled Hubble constant. 

Several solutions to the 5O(10) SUSY GUT dark matter problem have been proposed 
in Refs. [THE]- Here, we will concentrate on the most attractive one: that the dark matter 
particle is in fact not the neutralino, but the axino a. Axino dark matter occurs in models where 
the MSSM is extended via the Peccei-Quinn (PQ) solution to the strong CP problem|19j. The 
PQ solution introduces a spin-0 axion field into the model; if the model is supersymmetric, 
then a spin-| axino is also required. It has been shown that the a state can be an excellent 
candidate for cold dark matter in the universe [20]. In this paper, we will find that 5O(10) 
SUSY GUT models with an axino DM candidate can 1. yield the correct abundance of CDM 
in the universe, 2. avoid the gravitino/BBN problem and 3. have an compelling mechanism for 
generating the matter-antimatter asymmmetry of the universe via non-thermal leptogenesis. 

2 The gravitino problem 

An affliction common to all models with gravity mediated SUSY breaking (supergravity or 
SUGRA) models is known as the gravitino problem. In realistic SUGRA models (those that 
include the SM as their sub- weak-scale effective theory), SUGRA is broken in a hidden sector 
by the superHiggs mechanism, which induces a mass for the gravitino G, which is commonly 
taken to be of order the weak scale. The gravitino mass ttlq ends up setting the mass scale for 
all the soft breaking terms, so then all SSB terms end up also being of order the weak scale. 

The coupling of the gravitino to matter is strongly suppressed by the Planck mass, so the G 
in the mass range considered here (m^ ~ 5 — 20 TeV) is never in thermal equilibrium with the 
thermal bath in the early universe. Nonetheless, it does get produced by scatterings of particles 
that do partake of thermal equilibrium. Thermal production of gravitinos in the early universe 
has been calculated in Refs. [21] , where the abundance is found to depend naturally on m^, and 
on the re-heat temperature Tr at the end of inflation. Once produced, the Gs decay into all 
varieties of particle-sparticle pairs, but with a lifetime that can exceed ~ 1 sec, the time scale 
where Big Bang nucleosynthesis (BBN) begins. The energy injection from G decays is a threat 
to dis-associate the light element nuclei which are created in BBN. Thus, the long-lived Gs 
can destroy the successful predictions of the light element abundances as calculated by nuclear 
thermodynamics. 

The BBN constraints on gravitino production in the early universe have been calculated 
by several groups [22]. The recent results from Ref. [23] give an upper limit on the re- heat 
temperature as a function of tuq. The results depend on how long-lived the G is (at what stage 
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of BBN the energy is injected), and what its dominant decay modes are. Qualitatively, for 
rriQ ~ 5 TeV, T R ~ 10 6 GeV is required; if this is violated, then too many G are produced in 
the early universe, which detroy the 3 He, 6 Li and D abundance calculations. For txiq. ~ 5 — 50 
TeV, the re-heat upper bound is much less: T R ~5x 10 7 — 10 9 GeV (depending on the A He 
abundance) due to overproduction of 4 He arising from n «-> p conversions. For mQ ~ 50 TeV, 
there is an upper bound of Tr ~ 5 x 10 9 GeV due to overproduction of Z\ LSPs due to G 
decays. 

Solutions to the gravitino BBN problem then include: 1. having mg ~ 50 TeV but with 
an unstable Z\ (no Tr bound), 2. having a gravitino LSP so that G is stable or 3. keep the 
re-heat temperature below the BBN bounds. We will here adopt solution number 3. In the case 
of SO (10) SUSY GUT models, with ttiq ~ mi6 ~ 5 — 20 TeV, this means we need a re-heat 
temperature T R ~ 10 8 — 10 9 GeV. 



3 Non-thermal leptogenesis 

The data gleaned on neutrino masses during the past decade has lead credence to a particular 
mechanism of generating the baryon asymmetry of the universe known as leptogenesis |24j. 
Leptogenesis requires the presence of heavy gauge singlet Majorana right handed neutrino 
states ifjpfc(= Ni) with mass M^. (i = 1 — 3 is a generation index). The JVj states may be 
produced thermally in the early universe, or perhaps non-thermally, as suggested in Ref. [25] 
via inflaton <fi — > iVjiVj decay. The JVj may then decay asymmetrically to elements of the 
doublets- for instance Y(N\ — > h^e~) ^ T(N\ —> h~e + )- owing to the contribution of CP 
violating phases in the tree/loop decay interference terms. Focussing on just one species of 
heavy neutrino Aq, the asymmetry is calculated to be|26j 



r(JVi -> £+) - T(N 1 -> £-) 3 M A 



Ni 



T N 87T V, 
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m U3 S eff , (3) 



where m V3 is the heaviest active neutrino, v u is the up-Higgs vev and 5 e ff is an effective CP- 
violating phase factor which may be of order 1. The ultimate baryon asymmetry of the universe 
is proportional to e, so larger values of lead to a higher baryon asymmetry. 

To find the baryon asymmetry, one may first assume that the N\ is thermally produced in 
the early universe, and then solve the Boltzmann equations for the B — L asymmetry. The 
ultimate baryon asymmetry of the universe arises from the lepton asymmetry via sphaleron 
effects. The final answer [27], compared against the WMAP-measured result — ~ 0.9 x 10 -10 
for the baryon-to-entropy ratio, requires ~ 10 10 GeV, and thus a re-heat temperature 
T R ~ 10 10 GeV. This high a value of reheat temperature is in conflict with the upper bound 
on T R discussed in Sec. EJ In this way, it is found that generic SUGRA models are apparently 
in conflict with leptogenesis as a means to generate the baryon asymmetry of the universe. 

If one instead looks to non-thermal leptogenesis, then it is possible to have lower reheat 
temperatures, since the N\ may be generated via inflaton decay. The Boltzmann equations for 
the B — L asymmetry have been solved numerically in Ref. [28]. The B — L asymmetry is then 
converted to a baryon asymmetry via sphaleron effects as usual. The baryon-to-entropy ratio 
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is calculated in [22], where it is found 

where is the inflaton mass. Comparing calculation with data, a lower bound Tr ~ 10 6 GeV 
may be inferred for viable non-thermal leptogenesis via inflaton decay. 

4 Axino dark matter 

The sparticle mass spectrum described in Sec. [1] is characterized by 5 — 20 GeV scalars, but 
very light gauginos, with a \i parameter of order 1-2 TeV. As a consequence, the neutralino Z\ 
ends up being nearly pure bino. Since all the scalars are quite heavy, the predicted neutralino 
relic abundance ends up being very high: the calculation of Refs. [T8"j [T4"j find values in the 
range ^^Ji 2 ~ 10 2 — 10 4 , which is 3 — 4 orders of magnitude beyond the WMAP-measured 
abundance. 

A solution was advocated in Ref. [Hj that in fact the Z\ state is not the LSP, but instead 
the axino a makes up the CDM of the universe. The axino is the spin- 1/2 element of the axion 
supermultiplet which is needed to solve the strong CP problem in supersymmetric models. 
The axino is characterized by a mass in the range of keV-GeV. Its couplings are of sub-weak 
interaction strength, since they are suppressed by the Peccei-Quinn symmetry breaking scale 
f a , which itself has a viable mass range 10 10 — 10 12 GeV. While the axino interacts very feebly, 
it does interact more strongly than the gravitino. 

If the a is the lightest SUSY particle, then the Z-y will no longer be stable, and can decay 
via Z\ — > aj. The relic abundance of axinos from neutralino decay (non-thermal production, 
or NTP) is given simply by 

since in this case the axinos inherit the thermally produced neutralino number density. Notice 
that neutralino-to-axino decay offers a mechanism to shed large factors of relic density. For a 
case where mg- ~ 50 GeV and Q^h 2 ~ 1000, as can occur in 5*0(10) SUSY GUTs, an axino 
mass of less than 5 MeV reduces the DM abundance to below WMAP-measured levels. 

The lifetime for these decays has been calculated, and it is typically in the range of r(Z\ — > 
aj) ~ 0.03 secj20]. The photon energy injection from Zi — > S7 decay into the cosmic soup 
occurs well before BBN, thus avoiding the constraints that plague the case of a gravitino 
LSP[29j. The axino DM arising from neutralino decay is generally considered warm or even 
hot dark matter for cases with m« ~ 1 — 10 GeV[3T]. Thus, in our Yukawa-unified scenario, 
where mj ~ 80 MeV, we always get warm DM from neutralino decay. 

Even though they are not in thermal equilibrium, axinos can still be produced thermally in 
the early universe via scattering processes. The axino thermally produced (TP) relic abundance 
has been calculated in Ref. [2U1 ED] , and is given by 
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where g s is the strong coupling evaluated at Q — Tr and N is the model dependent color 
anomaly of the PQ symmetry, of order 1. The thermally produced axinos qualify as cold dark 
matter as long as m a ~ 0.1 Me V [201 130] - 

5 A consistent cosmology for axino DM from 50(10) 
SUSY GUTs 

At this point, we are able to check if we can implement a consistent cosmology for SO (10) SUSY 
GUTs with axino dark matter. Our first step is to select points from the 50(10) parameter 
space Eq. [T] that are very nearly Yukawa-unified. In Ref . [2] , Yukawa unified solutions were 
searched for by looking for R values as close to 1 as possible, where 

= max(f t , f b , f T ) 
min(f t , f b , f T ) 

where the f t , fb and f T Yukawa couplings were evaluated at M GUT . Thus, a solution with 
R = 1.05 gives Yukawa unification to 5%. 

We would like solutions where the axino DM is dominantly CDM. For definiteness, we will 
insist on Q^ TP h 2 ~ 0.01, while Qj p h 2 = 0.1. Thus, in step 1., we select models from the 
random scan of Ref. [H] that have R < 1.05, and m\Q : 5 — 20 TeV. In step 2., from the known 
value of and Q^^h 2 , we next calculate the axino mass needed to generate Q^ TP h 2 = 0.01 
according to Eq. [5J In step 3, we plug m„ into Eq. El where we also take g s = 0.915 (the running 
g s value at ~ 10 6 GeV), and PQ scale f a /N = 10 12 GeV. By insisting that VlT p h 2 = 0.1, we 
may calculate the value of Tr that is needed. 

Our results are plotted in the m a vs. Tr plane in Fig. [fl and occupy the upper band 
of solutions. In this plane, solutions with Tr ~ 3 x 10 7 — 5 x 10 8 GeV are allowed by the 
gravitino constraint (with ttIq ~ 5 — 20 TeV) and BBN. Solutions with Tr ~ 10 6 GeV can 
generate the matter-antimatter asymmetry correctly via non-thermal leptogenesis. Solutions 
with m a ~ 10~ 4 GeV give dominantly cold DM from TP of axinos. Solutions with mie > 15 
TeV are denoted by filled (turquoise) symbols, while solutions with mi 6 < 15 TeV have open 
(dark blue) symbols. 

We see that a variety of points fall in the allowed region. These points give rise to a 
consistent cosmology for 50(10) SUSY GUT models! Of course, there is some uncertainty in 
these results. We can take higher or lower values of the PQ breaking scale, higher or lower 
fractions of Q^ TP h 2 , and the Tr upper (and lower) bounds have some variability built into them. 
As an example, the lower band of solutions is obtained with Q^ TP h 2 = 0.03, Q^h 2 = 0.08 and 
fa/N = 5 x 10 11 GeV. In this case, some of the previously excluded solutions migrate into the 
allowed region to give a consistent cosmology with somewhat different parameters. 

6 Conclusion 

Our main conclusion can be summarized briefly. For Yukawa unified supersymmetric models, 
as expected in 5O(10) SUSY GUT models, we find one can implement a consistent cosmology 
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Figure 1: Plot of Yukawa unified solutions with R < 1.05 and 5 TeV < mie < 20 TeV in the ma vs.Tr 
plane. The upper band of solutions has VL^ TP h 2 = 0.01, ti~ p h 2 = 0.10 and f a /N = 10 12 GeV, while 
the lower band of solutions has n£ TP h 2 = 0.03, nj p h 2 = 0.08 and f a /N = 5 x 10 11 GeV. 

including the following: 1. BBN safe mass spectra owing to the multi-TeV value of rnie, which 
arises in SUGRA models from a multi-TeV tUq 2. a WMAP-allowed relic density of CDM that 
consists dominantly of thermally produced axinos, and 3. the re- heat temperature needed to 
fulfill the relic density falls above the lower bound required by non-thermal leptogenesis, and 
below the upper bound coming from gravitino/BBN constraints. 

We feel that the fact that Yukawa unified SO (10) SUSY GUT models pass these several 
cosmological tests makes them even more compelling than they were based on pure particle 
physics reasons. In any case, with a spectrum of light gluinos, charginos and neutralinos, 
they should easily be tested by experiments at the CERN LHC[32] even with low integrated 
luminosities of just ~ 0.1 fb^ 1 . 
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